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We investigate permutation-invariant continuous variable quantum states and their covariance matrices. We 
provide a complete characterization of the latter with respect to permutation-invariance, exchangeability and 
representing convex combinations of tensor power states. On the level of the respective density operators this 
leads to necessary criteria for all these properties which become necessary and sufficient for Gaussian states. 
For these we use the derived results to provide de Finetti-type theorems for various distance measures. 



I. INTRODUCTION 

A fc-partite density operator on (C d )® k is called n- 
exchangeable if it is the partial trace of a permutation- 
invariant state on (C d )® n . According to the quantum de 
Finetti theorem, such states can be approximated by convex 
combinations of fc-fold product states with error bounded by 
0(— ). This result has various applications in quantum infor- 
mation theory (e.g., in security proofs for cryptographic pro- 
tocols flO] or the justification of mean-field approaches ^) 
and relatives (the "monogamy" of entanglement (cf. ||4j,|5|]) or 
the appearance of a local classical description under symme- 
tries S, 0,1). 

Unfortunately, such de Finetti-type statements for finitely 
exchangeable states are no longer true when the Hilbert space 
C d of an individual ("local") system is replaced by an infinite- 
dimensional Hilbert space — as it happens for instance when 
dealing with fields of light. Indeed, as shown in f9[], there are 
?i-exchangeable states on (C d )® fe whose distance from any 
convex combination of product states is at least H(^). This 
outrules the possibility of there being a de Finetti theorem for 
all finitely exchangeable states on infinite-dimensional Hilbert 
spaces. 

It is, however, still possible to find de Finetti-type results for 
infinite-dimensional systems by considering restricted classes 
of states or observables J24ll . An example are the coherent cat- 
states for which a de Finetti theorem was obtained by D'Cruz, 
Osborne and Schack lHol] : an "exponential" version thereof is 
provided in lUTll . In the present work we will first restrict the 
set of observables and consider covariance matrices of canon- 
ical field operators (such as position and momentum or the 
quadratures of a field of light). In the second part we will 
then use these results in order to derive de Finetti-type theo- 
rems for a restricted class of states, namely Gaussian states. 
These play an important role in quantum optics, as coherent, 
squeezed and their thermal states are all Gaussian and remain 
so under the action of quadratic Hamiltonians and homodyne 
measurements. 

Symmetry and exchangeability for infinite-dimensional 
systems are defined analogously to the finite-dimensional 
case: We substitute the local Hilbert spaces C d by the Hilbert 
space Hd '■= Ti® d of d harmonic oscillators or modes with 
H = C 2 being infinite-dimensional. The symmetric group 
S n acts on the n-fold tensor product TL® n by permuting the 



factors. A state p on Tif n is called symmetric if it is in- 
variant under this action, i.e., if irpir^ = p for all n e S„. 
States tr w ®„-fc/9 on Ttf k obtained by tracing out n — k "lo- 
cal" systems of a symmetric state p on Ttf n are called n- 
exchangeable. Infinitely exchangeable states on Tt® k are 
those that are n-exchangeable for any n > fc; according to 
the infinite-exchangeability de Finetti theorem lfl2l [l3i [fill , 
which also holds on infinite-dimensional Hilbert spaces, these 
are all states of the form J T® k dm{r), i.e., convex combina- 
tions of fc-fold product states r® fe , where r are states on the 
single system Hd- We call such convex combinations power- 
states. 

We will provide some prerequisites in the next section and 
then give a detailed discussion of permutation-invariant co- 
variance matrices. This will yield necessary criteria for the 
respective density matrices which become sufficient as well 
when considering Gaussian states. De Finetti theorems for 
the latter are then provided in the second part. 



II. PRELIMINARIES 

This section will introduce the notation and recall basic re- 
sults (cf. lfl5ll ) on covariance matrices and the corresponding 
continuous variable states. Consider a system of n sites, each 
of which with d canonical degrees of freedom (modes). Such 
a system can be described using a symplectic vector space, the 
phase space, V = M. 2d <g> K" of dimension 2d ■ n. We choose 
canonical coordinates such that the symplectic matrix a has 
the form 



a = 



I,, 



Od 



(1) 



where 0^ and Id are the d x <i-zero and identity matrix re- 
spectively. Occasionally, we will denote with a slight abuse 
of notation the symplectic matrix for a subsystem (less than 
?i sites) by a as well. The block structure in Eq.fl]) reflects 
the grouping into position- and momentum coordinates. The 
corresponding quantum system with Hilbert space 7if n = 
(£ 2 )® dn is characterized by canonical position- and momen- 
tum operators {Rf}^Li for each site i = l,...,n; here 
{Rf, R d+a } for 1 < a < d are the position- and momentum- 
operator of the a-th mode of the i-th site, acting non-trivially 
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only on the i-th factor of the 77-fold tensor product space 7if n . 
These operators satisfy the canonical commutation relations 
(CCR) 

i[R%,R?] = a {a , kUM l. (2) 

Real linear transformations S : V — ► V which preserve 
the CCR (i.e., SaS T = a) are called symplectic. Affine 
transformations are obtained from the unitary Weyl operators 
W(£) := e 1 ^ aR rf25ll which give rise to displacements: 

W(£)B?W(tf = B? + &l. (3) 

These only alter the first moments of a state, collected in the 
displacement vector D := tr(pR), while keeping untouched 
the covariance matrix (CM) T with entries 

r(a,i),0Sj) = tr ({Rf - Df ■ 1, Bfj - Df ■ l}+p) . 

Since the Weyl operators act locally in the sense that they 
factorize with respect to local systems we can w.l.o.g. typi- 
cally assume D = 0. Positivity of the density matrices implies 
that any CM has to satisfy Y > ia. Conversely, every real 
symmetric matrix fulfilling this requirement is a valid CM. In 
particular, there always exists a Gaussian state which is, up to 
displacements, completely characterized by this CM. 

III. CHARACTERIZATION OF COVARIANCE MATRICES 

We are now ready to characterize covariance matrices 
with respect to permutation-invariance, exchangeability and 
power-state property. Clearly, we could express all results be- 
low in terms of Gaussian states. However, we want to logi- 
cally separate the treatment of CMs, since it is neither neces- 
sary for the state to be Gaussian nor does the corresponding 
density matrix have to share the full symmetry of the CM. 

We will say that a CM is a power-state CM if it corresponds 
to a mixture of states with product CMs of the form 7 <g> 1„. 
Similarly, we say that it is n-exchangeable if it is a sub-block 
of a permutation-invariant CM of n sites, and we call it sepa- 
rable if it corresponds to a mixture of states with product CMs 
oftheform0" =1 7j. 

The main ingredient of our analysis is the following one-to- 
one correspondence between permutation-invariant CMs and 
bipartite product CMs: 

Proposition 1 (Symmetric CMs). Let X be any real orthogo- 
naln x 77 matrix for which Xii = 1/ ^/n for alii — 1, ... ,77. 
Then S = \2d ® X is a symplectic transformation on V = 
M. 2d (g) M n which gives rise to a one-to-one mapping between 
permutation-invariant CMs YonV and pairs of CMs E,Fo« 
R 2d such that 

n 

r = r(E,F) = S(E0 0F)S T . (4) 

Proof. Evidently, S is a symplectic transformation and there- 
fore maps CMs onto CMs. In fact, it is a passive transforma- 
tion, i.e., number preserving. We first show that every pair of 



CMs E, F leads to a permutation-invariant Y. For this we have 
to show that the entries of the latter are site-independent in the 
sense that 

r (0,0,(73,0 = A*,/3 , ^( a ,i),(i3,j) = B a ,p for i ^ j . (5) 

The symmetric 2d x 2d matrices A and B encode the local 
(on-site) and non-local (inter-site) correlations, respectively. 
With some abuse of Dirac notation Eq.© can be written as 

Y = E®X\l){l\X T + F® {1- X\1)(1\X T ). (6) 

Exploiting the condition on the first column of X this leads to 

1 1 

F(o,0,(/3,.0 = ~^ a - fi + ~ n^"'' 3 ' ^ 

which is indeed permutation-invariant. This immediately im- 
plies that every permutation-invariant CM characterized by 
A, B can be obtained in the above way by choosing 

E = A + (77 - 1)B and F = A - B , (8) 

□ 

We will now use this representation in order to characterize 
power-state and exchangeable CMs: 

Proposition 2 (Power-state CMs). Let T(E, F) be a 
permutation- invariant CM on M. 2d <X> R". Then the following 
are equivalent: ( i) Y is a power-state CM, ( ii) Y is a separable 
CM, and (Hi) there exists a CM 70/d modes such that E > 7 
and F > 7. 

Proof. Clearly, (i) implies (ii). The remaining proof closely 
follows the derivation of the separability criterion for Gaus- 
sian states in lfl6ll . There it is shown that for every separa- 
ble CM T there exist CMs 7 l such that Y > 0" =1 Ji. Av- 
eraging over all permutations we obtain Y > 7 <g> 1„ with 
7 = i ^\ 7j. Applying the symplectic transformation of 
PropQ~fto this inequality leads to E, F > 7 (since the r.h.s. 
has E = F = 7) an thus (iii) is equivalent to (ii). Finally, 
again following Illfll . T > 7 ® 1„ implies that Y is the CM of 
a mixture of states with CM 7 ® l n . □ 

Proposition 3 (Exchangeable CMs). Let Y be a permutation- 
invariant CM on M. 2d <g) M fe , k < n. With the notation from 
Prop\l\Y is 77- exchangeable if and only if A + (77 — 1)B > ia 
or equivalently E=-^E+(1 — ^)F with E being any CM of 
d modes. 

This follows in a straightforward way from PropQ]and the 
relation Eq.®. It is remarkable that Prop.[3]gives a simple cri- 
terion for deciding 77-exchangeability for any CM. Note that 
such a criterion is not known for general exchangeable states 
on finite-dimensional Hilbert spaces (cf. Il7ll ). 

We also point out that according to Prop. [3] Y admits a 
permutation-invariant super-system of 77 sites for all 77 > k 
if and only if E > F; this is a stronger condition than the one 
for power-state CMs. 
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IV. DE FINETTI THEOREMS FOR GAUSSIAN STATES 

We now step back from phase space to Hilbert space and 
derive upper bounds on the distance between exchangeable 
Gaussian states and the set of power-states. Note that the con- 
ditions for a Gaussian state to be exchangeable or a power- 
state are precisely those given in the previous section. We 
begin with the single-mode case (d = 1) and recall that the 
trace norm (L 1 -norm) is defined by j| Aj|i = trV At A. 

Theorem 1 (Single-mode de Finetti). Let pr(E,F) be 
an n- exchangeable Gaussian density matrix describing a 
permutation-invariant k < n mode system. Define E' = 
^-^E. Then pr(E'.F) is a Gaussian power-state, and 



|Pr(E,F) - Pr(E',F) || x 



Proof. The matrix E' defines a valid CM, since multiplying by 
a factor A := — > 1 preserves the condition E > ia. Us- 
ing n-exchangeabiliry and the CM E appearing in Propf3]we 
can write E' = t^-jtE + F. Observe that E is a positive semi- 
definite operator, which is implied by E > ia (via complex 
conjugation and averaging). We conclude that E' > F. Ac- 
cording to Propf2] this implies that /Or(E'.F) is indeed a Gaus- 
sian power-state. 

To compute the distance between pr(E',F) an d Pr(E.F)> we 
use the tensor product form of the states in the basis of PropQ] 
together with the stability property || p®a — p'®a \ = \\p—p'\\ 
of the trace norm, so that 



|Pr(E,F) - Pr(E',F) 



PE - PE> 



(9) 



To evaluate this further we use that there is a canonical basis 
where E = sl2 for some s > 1 and thus E' = AsJ.2. The cor- 
responding density matrices are now simultaneously diagonal 
in Fock state basis (labeled by I = 0, ... , 00) with eigenval- 
ues pe(s) := 2(s — 1) /fa + l) e+1 and pi(Xs) respectively. 
Hence, we proceed with 12611 



\Pe - Pe'Hj 



J2\Ms)-MH (io) 

/ As - _ / a-i y+i' 
VAs + J Vs + 1/ 



max 2 



In order to obtain a bound which is independent of s we con- 
sider the worst case, i.e., we take the supremum over s > 1 
which is achieved for s = 1 1I27I1 for which 1 = attains the 
maximum so that finally 



Pe - pe 



1 A- 1 

I^ 2 aTT 



□ 



Let us now consider Gaussian systems with several modes 
per site (d > 1). In principle we could follow the same 
route as above. However, computing the trace-norm dis- 
tance directly appears to become cumbersome due to its non- 
additivity with respect to tensor products. For this reason 



we choose two different figures of merit for measuring the 
distance to the set of power-states: the relative entropy ifTill 
S(p,p') = tr[p(logp - log/?')] and the fidelity F(p,p') = 



try' y/pp'^p. Both have benign behavior under taking tensor 
products (they are additive and multiplicative, respectively) 
and both yield bounds for the trace norm distance since 

-\\P P\\i < I 2 -2F(p,p') 2 • 

Theorem 2 (Multi-mode de Finetti). Let Pr(E,F) be 
an n- exchangeable Gaussian density matrix describing a 
permutation-invariant system ofk<n sites of d modes each. 
Then the Gaussian power-state Pr(E'.F) with E' = ^3^E sat- 
isfies 



(n- ft/2 



S(/>T(E 1 F),0r(E' 1 F)J < log^ 



F (pr(E,F),Pr(E',F)) > 



n — k 
1 - fc/2 



d/2 



(ID 



(12) 



Proof. It is clear from the arguments in Thm[TJthat r(E', F) 
indeed corresponds to a Gaussian power-state. We again ex- 
ploit that the distance between pr(E,F) an d Pr(E',F) equals that 
between p£ an /?e' ■ Moreover, we can simultaneously diago- 
nalize the latter two states such that p^ = <3) a=1 p(s a ) where 
each p(s a ) = J2"eLo P(( s a)\ty (^| is a thermal Gaussian state 
with CM s a l2 and the s a > 1 are the symplectic eigenval- 
ues of E. For E' we only have to replace s a by A.s Q with 
A = n/(n — ft). For the relative entropy, summing up the 
series leads to 



£(/>r(E,F),/>r(E',F)) = y^S(p(s a ),p(Xs a )) (13) 

a=l 

1 d 

a=l 

+ (s a + 1) log(As Q + 1) - (s a - 1) log(As a - 1)) . 

For every A > 1 each summand is a decreasing function in s a 
such that the supremum is again obtained for s a = 1, i.e., 

sup S(p(s),p(Xs)) =logii^ 

s>l z 

which concludes the proof for the relative entropy upon in- 
serting A. In a similar vein we can evaluate the fidelity: 

d 

Pr(E,F),Pr(E',F) ) = " lF(p(s a ),p(Xs a )), (14) 

Q = l 



F{p(s),p(^s)) = ^2^p e (s)pe(\s) 

1=0 

2 



V(s + 1)(A S + 1)- V(s-1)(A S -1) 



(15) 



Again the worst case, now the infimum over s > 1, is attained 
at s = 1 which leads to the desired result. □ 
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V. CONCLUSIONS & OPEN PROBLEMS 

The bounds given in ThmQ] and Thm|2] are tight for the 
given class of ansatz states and the proofs allow for any ex- 
plicitly given CM E to compute the exact distance to this class. 
Whether the chosen class is optimal or yields at least an op- 
timal scaling remains open. Similarly, the question for which 
figure of merit a power-state closest to a given n -exchangeable 
Gaussian state can again be chosen to be Gaussian remains 
open as standard arguments, e.g. based on the central limit 
theorem (CLT) 11911 . do not immediately apply. Some kind of 
marriage between the CLT and the de Finetti theorem would 
also be desirable in order to extend the latter in a reasonable 
way beyond the class of Gaussian states. 

Another interesting direction of future research could be to 
use more general approximating states than Gaussian power- 



states in the de Finetti theorem. In this way it might be pos- 
sible to obtain an exponentially small error as in the "almost- 
product" de Finetti theorem for finite-dimensional systems (H 
(see also |[Tlll20ll ). This would be important for applications 
to continuous-variable quantum key distribution, for exam- 
ple & 
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